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SUMMARY 


The  method  of  "caustics"  is  investigated  and  analyzed  as  a  means  of 
evaluating  stress  intensity  factors,  Kj,  in  cracks  radially  emanating 
from  holes  drilled  in  aluminum  and  steel  specimens.  The  results  from 
cracks  in  open  holes  and  in  specimens  under  uniaxial  tension  loads, 
are  compared  with  results  obtained  analytically  (as  described  by 
Paris  and  Sih  in  ASTM  381-Ref  5).  Effects  of  compressive  residual 
stresses  in  reducing  are  demonstrated. 

The  classical  analyses  of  stress  concentrations  in  the  vicinity  of 
a  crack  tip  have  led  to  the  determination  of  stress  intensity  factors 
that  are  useful  in  predicting  design  allowable  stresses  in  riveted 
joints.  These  analytical  predictions  have  now  been  related  to  direct 
measurement  of  surface  deflections  by  the  method  of  caustics,  which 
employs  laser  reflection  under  controlled  conditions.  The  caustic  is 
the  trace  of  the  singular  reflection  zone  in  the  image  plane  of  an  ex¬ 
perimental  apparatus,  and  its  dimensions  are  a  very  sensitive  measure 
of  surface  slopes.  This  method,  originally  devised  by  P.S.  Theocaris 
in  1972,  has  been  successfully  employed  in  a  variety  of  applications. 

In  this  report,  we  present  initial  results  of  its  application  to  the 
investigation  of  crack  propagation  under  load  as  it  is  affected  by  re¬ 
sidual  compressive  stress  distributions  created  by  dynamically  installed 
interference  rivets.  When  compared  with  analytical  prediction  and  the 
results  of  fatigue  tests,  we  hope  that  eventually  it  will  be  possible  to 
define  in  advance  operating  domains  within  which  crack  growth  will  not 
occur. 

Several  clear  advantages  were  offered  by  the  method  of  caustics. 

For  example,  classical  analytical  models  imply  that  Kj  =  0  when  the  load 
is  zero.  Yet  we  know  that  in  any  crack  the  residual  stress  field  at 
the  tip  will  always  mean  that  Kj  /  0.  Since  the  method  of  caustics 


shows  large  changes  in  light  reflection  from  small  surface  deformations 
near  the  tip  of  the  crack,  it  measures  Kj  accurately  enough  to  determine 
Ki  at  no  load.  Another  advantage  the  caustic  gives  is  the  determination 
of  the  load  at  which  plastic  deformation  begins.  Again,  current 
analytical  models  cannot  predict  such  loads,  but  the  sudden  enlargement 
of  the  caustic  is  a  very  good  indication  of  incipient  plastic  deformation. 

Finally,  the  holes  in  aluminum  and  steel  specimens  were  cold  worked 
using  a  dynamically  expanding  rivet.  We  observed  that  whenever  the 
radius  at  which  the  elastic-plastic  boundary  was  expanded  beyond  the 
crack,  the  caustic  disappeared;  but,  if  the  elastic-plastic  radius  was 
smaller  than  the  crack,  the  caustic  increased  in  size.  We  believe  that 
the  method  of  caustics  can  be  used  with  great  advantage  to  evaluate  the 
effectiveness  of  cold-worked  holes  in  riveted  joints. 
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damage  tolerance  criteria  have  been  used  in  recent  wars  to  evaluate 
the  economic  lift'  oh  uifvi'Jht  structures.  lhe  current  Mil -A-S3444 
damage  Tolerance  Spec i f icat ion  requires  such  criteria  to  be  used  in 
design.  They  are  based  on  the  premise  that  the  structures  can  with¬ 
stand  the  existence  oh  small  cracks  without  catastrophic  failures, 
tor  tiie  duration  oh  their  predicted  fatigue  life.  There  are  many 
critical  parts  in  aircraft  structures  for  which  premature  failure  can 
cause  the  withdrawal  of  the  aircraft  from  service  because  the  cost  of 
its  repair  will  exceed  the  cost  of  a  new  aircraft.  This  limit,  called 
the  economic  life  of  the  aircraft,  is  one  of  the  design  criteria 
used  to  minimice  total  life  cycle  costs.  One  of  the  most  common 
sources  of  such  failures  in  aircraft  is  in  the  riveted  joints.  Recently, 
i Ret  11,  we  showed  that  considerable  fatigue  life  enhancement,  as  well 
as  good  economic  life  and  lower  initial  cost  ,  can  be  achieved  if  the 
joints  are  riveted  with  residual  compressive  stresses.  Any  evaluation 
of  riveting  methods,  however,  that  provides  these  advantages,  involves 
the  evaluation  of  the  stress  intensity  factors,  and  the  residual  stress 
field  from  which  wo  can  evaluate  the  remaining  fatigue  life  of  damaged 
^cracked)  structures.  We  have  used  the  new  riveting  method  called 
Stress  Wave  Riveter  vSWRl  ^Ref  .')  to  achieve  a  considerable  degree  of 
tgtioue  life  enchancement  in  such  damaged  joints. 

An  experimental  method  of  evaluating  stress  intensity  factors  (kj) 
is  necessary  to  provide  us  with  the  means  to  evaluate  the  degree  of 
enhancement  due  to  a  particular  residual  stress  distribution  and  flaw 
sice  quickly  and  accurately  i refer  to  Appendix  Al.  We  selected  the  method 
of  caustics  to  explore  the  feasibility  of  measuring  kj .  Theocaris  (Refs  3 
and  4'  has  done  considerable  work  on  the  theory  and  applications  of  caustics. 
A  detailed  analysis  of  caustics  is  given  in  Section  2.  To  illustrate  how 
the  method  of  caustics  works,  we  analyzed  a  one  dimensional  model  of  a  crack 
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emanating  from  a  hole,  and,  by  using  the  principles  of  reflecting  rays,  we 
were  able  to  show  how  the  coalescing  of  light  rays  form  caustics. 

We  are  using  a  laser  beam  as  an  intense  collimated  monochromatic 
light  source.  From  an  optical  point  of  view  unfinished  metal  surfaces 
are  not  good  reflectors.  To  produce  good  caustics,  therefore,  we  must 
make  every  possible  effort  to  obtain  smooth  surfaces,  especially  by 
polishing  the  specimens.  Caustics  can  be  used  in  many  appl ications . 

We  are  interested  here  only  in  their  application  to  cracks.  In  par¬ 
ticular,  we  are  using  caustics  as  a  non-destructive  means  of  evaluating 
riveted  aluminum  and  steel  structures.  In  our  efforts  to  improve  fa¬ 
tigue  life,  and  to  reduce  the  cost  of  manufacturing  and  maintaining 
aircraft,  we  feel  that  diagnostics  such  as  the  method  of  caustics  to 
detect  and  evaluate  cracks  will  eventually  prove  to  be  practical. 

Our  work  consists  of  making  fatigue  specimens  with  holes.  Using 
cyclic  loads,  we  precrack  each  hole  to  a  specified  length.  The  speci¬ 
mens  are  polished  prior  to  cracking  in  accordance  with  standard  metal¬ 
lurgical  procedures.  The  caustics  for  each  cracked  hole  are  then 
measured  and  the  Kj  factors  evaluated  at  various  tension  loads. 

Using  the  SWR  method,  rivets  are  then  installed  with  various  degrees 
of  compressive  residual  stresses  (cold  work).  The  Kj  values  are 
measured  again  at  various  tension  loads.  Details  of  the  laboratory 
set-up,  test  results  and  evaluation  are  given  in  the  following 
sections  of  this  report.  The  results  show  clearly  that  there  are 
regimes  of  crack  length,  load  and  residual  stress  levels  in  which 
stress  intensity  factors  are  reduced,  and  in  some  cases,  eliminated 
completely. 
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2.  THEORETICAL  BACKGROUND 

THEORETICAL  EXPLANATION  OF  THE  METHOD  OF  CAUSTICS  FOR  MEASURING  K 

The  derivations  of  this  section  follow  the  methods  described  by 
P.S.  Theocaris  and  his  colleagues  (Ref  3).  We  have,  however,  developed 
the  uniaxial  case  and  derived  the  "ray  reflection"  method  to  elucidate  the 
particular  case  of  caustic  application  to  cracks. 

The  surface  deformation  around  the  tip  of  the  crack  is  given  by 
Eq.  (A-4)  from  the  Appendix 


*{t)=  2(2*£>1/2  U) 

where  K*  is  the  complex  stress  intensity  factor  and  t  =  x  +  iy. 

Considering  a  cross  section,  in  the  direction  of  the  crack,  with 
only  one  mode  of  crack  opening,  Eq  (1)  is  reduced  to 

2(2nx)l/2  ^ 

The  surface  deformation,  on  the  other  hand,  is  given  by 

Z  =  At  =  —  —  (<7t  +  a2)  (3) 


where  v  =  Poisson's  ratio 
E  =  Young's  modulus 

t  =  thickness 

and  oj  &  u2  are  the  principal  stresses  around  the  crack.  Paris  and  Sih, 
(Ref  5),  however,  related  Eq  (2)  &  (3)  by  the  following: 


Therefore 


o  ( +cr2  =44>(Z) 


Z  = 


ktK, 

E(27rx)1/5 


(4) 


Differentiating  (6)  we  obtain 


dz  _  1_  C 
dx  2  x3^2 


d2Z  _ _ 3  C 
dx2  4  x5/2 

The  geometrical  representation  of  caustics  is  shown  in  Fig.  2-1 


RB 1-0244-042 W _ _ _ 

Fig.  2-1  Reflection  of  Incoming  Ray  from  Surface  Adjacent  to  Crack 


The  angle  of  the  reflected  ray  is  given  by 


a  =  90  +  26  +  <p  . 

„  .  dZ  .  C 

6  =tan  Tt  =  tan  2^T 


=  tan" 


Thus, 


a  =  90  +  2  tan"1 


2?7r 


To  illustrate  the  meaning  of  Eq  (10),  let  us  consider  a  case  of  an  alumi¬ 
num  specimen,  0.125  in.  thick,  with  Kj  =  5  x  10^  psi/TrT. 

Hence  v  =  0.33 


E  =  107  psi 
t  =  0.125  in. 

Substituting  in  Eq  (5)  we  obtain 

C  =  0.8228  x  10'5 

Therefore,  using  Eq  (10)  we  find  the  following: 


0.011 

0.013 

0.014 

0.01435 

0.0145 

0.015 

0.016 

0.017 


90.7238 

90.6905 

90.6856 

90.6853 

90.685 

90.688 

90.691 

90.700 

90.740 


From  a  close  look  at  the  change  of  a  with  x,  we  see  that  between 
x  =  0.014  and  x  =  0.015,  Aa  changes  from  negative  to  positive.  In  fact, 


the  point  where  this  change  takes  place  is  located  at 

3  6  =  (p. 


Since  the  angles  0  and  <j>  are  very  small,  we  have 


from  which  we  obtain 


3  C  _  1_ 

2  x5/2  Z, 

Comparing  Eq  (8)  with  Eq  (11)  we  see  that  Act  changes  sign  at 

d2Z  _  _  1 
dx2  _  2Z( 


(ID 


(12) 


Solving  Eq  (11)  for  x 

(3  \ 2/6 

—  CZj )  .  (13) 

where  x..  is  the  innermost  point  inside  of  which  reflected  rays  diverge. 
For  the  numerical  case  we  considered  above  we  have 

Xj  !=0.  01435  in. 

or,  an  angle  0  =  0.1371°  and  a  =  90.6853. 

If  we  now  place  a  screen  at  a  distance  Zq  from  the  specimen  we 
have  the  geometry  shown  in  Fig.  2-2. 

Considering  the  triangle  ABC,  in  Fig.  2-2,  we  have 

W  =  x-(ZO  +  Z)tan(20  +  tp)  (13) 


or 


w  =  x  —  (Zq  +  Z) 


tan20  +  tan  (fi 
1  -  tan26  tan  <p 


Since  e  and <p  are  small  angles,  we  set  tan2e  =  2tane  =  2(dZ/dx)  and 
tan2o  tan<A  s  2tane  tan s  0;  tan^?=  x/Z^ . 

From  our  previous  discussion  we  know  that  below  the  value  (for 
example  x^  =  0.01435  in.)  reflected  rays  will  diverge.  These  rays  do 
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not  contribute  to  the  caustic  formation.  The  value  of  is  given  by 
Eq  (13).  For  x  >  x^ ,  however,  reflected  rays  converge.  The  condition 
of  convergence  of  two  or  more  rays, at  the  O' x'  axis  on  the  screen,  is  ex¬ 
pressed  by  setting  the  Jacobian  equal  to  zero.  In  our  one  dimensional 
model  this  condition  is  given  by 


Differentiating  Eq  (13)  we  obtain: 


dw 

dx 


1  -  2Z„ 


d2Z 

dx2 


Z A 
Z» 


=  0  . 
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Hence 


1  +  ^  =  -  2Zn 


d2Z 
d x5" 


and 


(14) 


^  _  1  -  _  2Z  -^2- 

z,  1  2Z°^?-  • 

Zo  +  Zi 

The  ratio  “  =  A  represents  the  amplification  factor  of  a 


Zi  tn 

converging  laser  beam  (Note  Z^  <  0).  Hence, 


2Z  o  d2Z 

=  ~  1 


\n  dx 

Substituting  from  Eq  (8)  we  have 


3  C 


4  Xj/2 


2Z0 


and 


5/2  _  3CZ^ 
1  2X„ 


(15) 


and 


x'=  r0  = 


|  3wtZnKt 
|  2E(2tt)1' 2 


2/5 


(16) 


Thus,  the  point  where  two  converging  rays  meet  in  the  O’x'  axes  is  given 
by 


x'  = 


Rearranging  Eq  (13)  we  have: 


w  =  x-2Z0  ^  -Z0tan<p 


and 


or 


W  =  X-Z0^7T 


ZflX 

z, 


=  x 


w  =  x 


Zt 

Znx 

Zt 


(17) 


(18) 

(19) 
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Substituting  from  Eq  (15)  we  have 


or 

and  finally,  w  =  5/3(Amx^)  (i.e.,  maximum)  and  w  =  1/3  Afflx.  (i.e. 
mum.)  Thus,  from  Eq  (15)  we  have: 

2E(2tt)1/2  ^ . . 


2E(27r),/2  f  D\in 


6 


MZ0\n  3vtZ0 

where  D  is  the  diameter  of  the  caustic  perpendicular  to  the  crack  axis 
and  D  =  6r0\„.  The  caustic  formation  can  be  shown  graphically  (see 
Fig.  2-3). 

To  evaluate  D,  it  is  necessary  to  carry  out  the  two  dimensional 
case  (Ref  3).  For  this  purpose,  consider  a  specimen  with  a  hole  from 
which  a  crack  is  growing  (see  sketch  below). 
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At  a  point  P  the  deflection  due  to  the  crack  is  given  by 


K  ^  c 

=  2(27r£  )*/ 2  =  2(2I7rr),/5 


K,  e' 


•i(oj/2) 


(23) 


where 


£  =  x  +  iy  =  relw 

The  half  distance  w  of  the  caustic  from  the  screen,  placed  Zg 
distance  away  is  given  by:  w  =  ?  -  (Zg  +  Z)  tan(2e  +  </>  ).  This  equa¬ 
tion  can  be  further  reduced  to, 

w=£  -2Z0tan 6  -Z0tan <p; 

using  Eqs  (3)  &  (4)  we  obtain 


w=£+Z0~  grad  [  4  Re  $(£)!- 

t  Lt  ^ 


(25) 


From  Eq  (13)  we  derive  the  Rei>U)  to  be, 

CO 

Kx  cos  — 

u(r>  wl  = 

The  grad  u  =  (au/ar)  +  i  ( 1/r )  Ou/d0),  since  s>U)  is  a  harmonic 
function.  Hence, 


(26) 


Kx 


Substituting  Eq  (27)  into  Eq  ( 25 ;  we  obtain 


,  U> 

.  u> 

7  u)  .  .  u>  \ 

(  cos- 

ismF) 

Kjlcosg-  +  ism-7 

("  2r3/i 

2ri/2  J  ' 

4(2tt),/2  vi/2 

(27) 


w  =  x'  +  iy'  =  —  r  elu 


/Zn  — Z,  \  .  c 

E(2AW! 


ZnKKr  e~Hw/2) 


or 


where 


X-  My-  {  - 1  +  |  (28) 


Considering  now  the  (x,  y)  coordinates  where  t  =  x  +  iy,  and 
(x‘,  y')  coordinates  of  the  screen,  we  know  that  there  is  not  one  to 
one  correspondence,  since  distinct  points  on  the  Z(x,  y)  converge  into 
a  single  point  in  the  (x1,  y')  plane.  This  statement  is  expressed  by 

11 


the  Jacobian  J  =  0 
or 


where. 


9(o.'t,  o-'v) 

8(x,  y) 


9ta„ 

9ia„ 

9x 

9y 

dil)y 

9U)y 

9x 

9y 

-  2Z0 


wy  Xmy  2Z0  ) 

and  f(x,  y)  =  -  (4vt/2E)Re-i’(£)  represents  the  surface  deformation  due  to 
the  crack  formation. 

Differentiating  each  of  the  Eq  (30),  substituting  into  (29)  and 
evaluating  the  determinant  we  obtain  the  following: 

‘0 

°r  v(0^0)-zi[(0-)(0)-(e^)!]=o  (31> 


*  -  «.».  ^ ^ - (sfij  J  - O 

The  function  f,  however,  can  be  considered  harmonic  (Ref  5)  and 
therefore  the  Cauchy-Rieman  relations  apply. 

Thus,  if 

f(£)  =  u(x,  y) +iv(x,  y),  where  J  =x  +  iy,  j  ( 


we  have 


9u  9v  ,  9v  9u 

9x  9y  an  9x  9y 


As  suggested  by  Paris  and  Sih,  however,  we  need  only  consider  the 
real  part  of  f(c). 

From  Eq  (32)  we  can  easily  prove  that  (a2u/ax2)  +  (32u/ay2)  =  0, 


and  since 


32f  _  32u  .  32v 
9£"  9x^  1  9x^ 
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and 


Furthermore  g»  .  .  0  .  hence  (§|) ‘  -  -  (g)  (gj) 

/ 92vV  /  &2U  y  (  (34) 

\  9x9y  / 

Substituting  Eq  (33)  and  (34)  into  Eq  (31)  and 

Mwr-- 


(35) 


From  Eq  (23)  we  obtain 


3 2  4>  =  3  KTr5^2 
dp~  8  (2ir),/2 

,  2 


and 


d2<fi  |  (3  K,r's/2  1 

8F|  =  l  8  l2^"| 


Substituting  into  Eq  (35)  we  have: 

r 5 / 2  _  3 _ K>vtZ0 

0  2  (27r)1^EXro  ' 

Equation  (37)  is  valid  only  at  the  point  where  Eq  (35)  holds  true. 
Substituting  in  Eq  (28)  we  have. 


(36) 


(37) 


w 


=  r0  Am  e‘"^-l+|e'‘(3“/2^ 


(38) 


The  value  of  rQ  is  equivalent  to  the  value  of  of  Eq  (13)  in 
the  uniaxial  case  and  it  represents  the  limit  that  separates  the  rays 
that  converge  from  the  rays  that  diverge. 


Thus, 


and 


:  r0  \  COSU)  +  -  COS 


y'  =  ro\i,  <  sinw  + 


3cu  ) 

2) 

2  .  3u>  ) 

3  Sin~ ] 


(39) 


At  the  point  where  x1  is  zero  we  have  y1  =  D/2,  while D  =  6rg  along 
the  y1  axis,  x'  =  0  is  equivalent  to 


2  3u>  n 

COS  u;  +  -  cos  — —  =  0  . 

O  c* 
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This  equation  is  satisfied  in  the  first  quadrant  at  to  =  75°.  Evaluat¬ 
ing  y'  at  to  =  75°  we  obtain  y1  =  1.582  r^.  Therefore,  D  =  3.164  rQXm 
from  which  we  have  5  =  3.164. 


Rearranging  Eq  (22)  we  obtain: 


K, 


1.671 E 


^tz0xm 


J7T 


(40) 


By  measuring  D  on  the  screen  at  a  distance  Zq,  we  obtain  the  value  of 

Kr 

When  we  analyzed  the  uniaxial  case  at  y  =  0,  we  found  the  interest¬ 
ing  relationship  expressed  by  Eq  (11)  independently  of  the  Jacobian,  by 
observing  the  values  of  the  angle  a.  Can  Eq  (11)  apply  at  any  angle, 
other  than  to  =  0?  Using  r3  instead  of  x  in  Eq  (11)  we  can  easily  find  the 
unique  solution  of  to  =  0°.  Hence,  3b  =  <j>  is  true  at  w  =  0  only. 

AN  ILLUSTRATIVE  EXAMPLE  OF  CAUSTICS 

Let  us  consider  a  set  of  coordinates  (Z,  x)  in  a  plane  (see  sketch 
below). 


A  crack,  with  its  center  at  x  =  a  can  be  represented  by  the  following 
equation . 


K2 

x  =  a±^j  ,  where  a > 0  and  x > 0 


For  a  crack  with  its  center  at  x  =  -  a  we  have 

K2 

x  =  -  a  *  — jy  ,  where  a  <  0  and  x  <  0  . 
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In  both  cases  K  is  a  constant  to  be  determined  experimentally. 


Let  us  now  consider  the  case  of  a  =  -  0.630  in.  For  the  branch 
x  <  a,  we  have. 


z 


K 

(a  -  x)1^2 


(41) 


and  for  the  branch  x  >  a 

K 

Z  "  "  (x“'aF* 


(42) 


Of  course,  at  x  =  a,  we  have  a  discontinuity,  a  characteristic  of  any 
crack. 

Differentiating  Eq  (41)  and  (42)  we  obtain 


and 


dZ  _  K  1 

dx  2  (a-x)3y/i  ’  X  a 


dZ 

dx 


K  1 
2  (x  -  a 


x>  a  . 


(43) 

(44) 


Let  a  converging  light  beam  with  a  focus  at  Z-j  =  -  1.875  in.,  be  di¬ 
rected  on  the  surface  that  contains  the  crack  along  (AB)  (see  sketch  below). 


-i  |  _  2^-?)  -2  |  +  tan'1  (  7-)  X  where  Z  «  Zi ; 


From  the  geometry  of  the  diagram  we  easily  deduce  that 
a?'  =  20  —  <p  =  -  2  tan" 

but  a  =  180  -  u*  (w‘  is  measured  clockwise,  and  w  counterclockwise) 
Hence,  u>  =  180  - 2  tan"1 1  —  ~  )  +  tan"1  1  ^  I  (47) 


when 


u 


=  80°,  zi  =  —l.  875  in. ,  a  =  -  0. 630  in. ,  and  x  =  -0.50  in. 


(These  values  were  measured  from  a  geometric  construction  of  caustics) 
Hence , 

1  875 

u,' =  180-0.)  =  100%  C^tan"1  -r— — -  =  75° 

U,  DU 


and 

Therefore 

or 


c o'  +  <p 

6  = - — ^  =  87.  5“  . 

Z 


tan#  - 


2(x  -  a):l/2 
K 


_  2(x-a):i/2  _  2(—  0.  50  +  0.  630)1*5  _  Q  Qn 


tan0 


tan  87° 


and  K  =  4.  91  x  10'3 

Thus,  the  angle  w  for  various  values  of  x  is  given  by  Eq  (47)  with 
K  =  4.91  x  10'3 


OJ 


or 

for 

and 


=  180  —  2  tan"1  {203.  04  (x  +  0. 630)*'b}  +  tan"1  ^ 


x> -0.630 


cu 


=  -  180+2tan"'{203.04(-(0.  63+x))K5}  +  tan"1^--^-5^ 


(48) 


(49) 


for  x< -0.630 

The  results  are  plotted  in  Fig.  2-5. 


X 
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Fig.  2-5  Formation  of  Caustics 
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A  caustic  obtained  from  a  precracked  open  hole  on  an  aluminum  speci¬ 
men  is  shown  in  Fig.  2-6.  In  this  example,  the  crack  is  long  and  we  can 
distinguish  the  straight  portion  that  is  already  separated  from  the  cir¬ 
cular  portion  which  represents  the  tip  of  the  crack.  In  this  example, 
the  crack  is  at  the  Z  axis  and 


It  represents  the  uniaxial  case  discussed  earlier  in  this  section.  Infor¬ 
mation  pertaining  to  the  specimen  is  given  in  Fig.  2-7.  The  value  of  D 
was  measured  and  Kj  was  evaluated  from  Eq  (40)  for  various  tension  loads 
o.  The  details  are  given  in  Section  3.  The  theoretical  estimate  of  Kj 
is  given  by  Kj  =  kj/rra  where 

a  =  crack  length 

k  =  depends  on  the  ratio  (a/d),  (d  =  hole  diameter) 
o  =  uniaxial  stress  on  specimen 

This  equation  was  taken  from  ASTM  381  by  Paris  and  Sih  (Ref  5). 

The  caustics  of  the  same  crack  under  various  loads  are  shown  in 
Fig.  2-8.  The  increase  of  the  large  circular  area  at  the  tip  of  the 
crack  with  load  corresponds  to  the  graph  shown  in  Fig.  2-7. 

CRACKS  IN  HOLES  WITH  RESIDUAL  STRESSES 

From  work  done  previously  (Ref  1)  we  Know  that  fatigue  life  in 
riveted  joints  increases,  if  the  fastener  generates  residual  stresses  in 
the  radial  and  circumferential  directions.  Some  fasteners,  like  Hi-Tique 
and  Taper-Lok,  generate  radial  compressive,  but  circumferentially  tensile 
stresses.  This  is  accomplished  through  radial  interference  by  installing 
a  pin  with  larger  diameter  than  the  hole.  However,  greater  increases  in 
fatigue  life  can  be  obtained  when  a  radially  expanding  rivet  plastically 
deforms  the  material  around  it,  so  that  upon  relaxation  of  the  installa¬ 
tion  forces,  compressive  residual  stresses  are  generated  in  both  the 
radial  and  circumferential  directions.  This  method  of  rivet  installation 
is  also  known  as  "cold  working  the  hole". 

The  degree  of  fatigue  life  enhancement  depends  on  the  extent  we  cold 
work  the  hole,  more  specifically,  how  large  the  radius  of  the  elastic- 
plastic  interface  becomes  during  plastic  deformation.  The  elastic-plas¬ 
tic  interface  is  defined  as  the  radial  position  (rp),  inside  which  we 
have  a  compressive  biaxial  field  and  beyond  which  the  circumferential 
stresses  become  tensile  (see  Fig.  2-9).  If  a  crack  emanating  from  the 
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STRESS  INTENSITY  FACTOR  -  KSt  v/lN. 


Fig.  2-7  Strew  Intensity  Factors  (ksi  >/Tn.).  Caustics  vs  Theory,  Specimen  3C 
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hole  exists  prior  to  the  cold  working,  with  length  a  <  rp,  a  most  desir¬ 
able  situation  exists  for  high  fatigue  life  enhancement,  since  the  bi¬ 
axial  compression  field  will  retard  the  crack  growth.  In  this  region, 
therefore,  we  expect  to  see  a  reduction  in  Kj  values,  in  comparison  with 
Kj.  values  in  precracked  open  holes,  until  the  applied  load  on  the  speci¬ 
men  exceeds  the  residual  stresses.  If  the  crack  length  a  >  r  ,  then  the 
initial  kj  values  will  be  greater,  after  the  cold  working  of  the  hole, 
than  the  Kj  values  obtained  with  the  precracked  open  hole,  but  will  re¬ 
main  constant  until  the  applied  loads  in  the  specimen  exceed  the  residual 
stresses. 

A  rivet  was  installed  in  the  precracked  hole  which  generated  com¬ 
pressive  residual  stresses,  not  extensive  enough  so  that  a  >  r  .  Fig¬ 
ure  2-10  represents  the  caustics  of  a  precracked  hole  (described  in  the 
previous  section)  under  load,  before  riveting  (open  hole).  Figure  2-11 
presents  the  caustic  after  cold  working.  Equation  (40),  which  was  de¬ 
rived  earlier,  to  relate  with  the  major  diameter  of  caustics  applies 
to  the  results  shown  in  Fig.  2-10,  but  it  does  not  apply  to  the  results 
shown  in  Fig.  2-11,  because  the  conditions  after  riveting  are  no  longer 
elastic.  Plastic  deformation  has  distorted  the  caustic  in  Fig.  2-11,  but 
it  did  not  eliminate  it.  In  other  cases,  where  the  cold  working  is  more 
extensive,  and  the  cracks  smaller  (i.e.,  a  <  rp)  the  caustic  has  com¬ 
pletely  vanished  after  riveting. 

Finally,  the  rivet  was  drilled  out.  The  residual  tension  hoop 
stress  is  now  larger,  causing  an  increase  in  the  size  of  the  caustic. 
Figure  2-12  presents  the  caustics  at  various  loads  for  a  precracked, 
cold  worked  hole  under  various  loads.  We  can  easily  see  that  the  sizes 
of  the  caustics  in  Fig.  2-12  under  corresponding  loads  with  the  open 
hole  results  of  Fig.  2-10,  are  larger  and  fatigue  results  can  easily 
prove  that  in  this  case  (i.e.,  where  a  >  r  )  cold  working  has  worsened 
the  initial  cracked  hole  conditions  for  fatigue  life. 

It  is  this  qualitative  comparison  of  experimental  results  in  cold 
worked  holes,  that  can  be  obtained  by  the  method  of  caustics  in  a  unique 
fashion.  Cold  working  holes  through  riveting  has  been  proven  effective 
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in  fa'iqur  I  re  ■  nh.ir.  .emont  (  i  .e.  ,  increar  .d  by  a  factor  of  20)  of 
riveted  joints,  but  care  must  be  taken  to  provide  the  degree  of  cold 
work  appropriate  for  the  crack  sizes  we  are  trying  to  retard.  The 
method  of  caustics  gives  us  a  means  to  evaluate  difficult  problems  that 
usually  are  not  easily  amenable  to  analytical  solutions.  In  addition, 
we  can  save  valuable  time  in  evaluating  conditions  before  repairs,  and 
provide  a  non-destructive  means  of  evaluation  after  the  repair  of  riveted 
joints. 


6.0  KSI  10.0  KSI 


3.  EXPERIMENTAL  SET-UP 


The  experimental  set-up  used  for  this  series  of  tests  of  the 
caustic  phenomenon  was  assembled  as  a  prototype  system  from  equipment 
available  in  Grumman  Aerospace  Corporation's  Experimental  Mechanics 
Laboratory.  This  set-up  is  shown  in  Figs.  3-1  and  3-2. 

The  set-up  consists  of  three  basic  components:  a  light  conditioning 
assembly,  a  specimen  loading  frame,  and  a  screen  and  camera  to  image 
and  record  the  caustic  display. 

The  light  conditioning  assembly  is  composed  of  a  platform  holding 
the  optical  components  mounted  on  a  combination  of  stages  which  allow 
controlled  motion  in  the  x-y-z  directions,  and  rotation  about  the  z 
axis.  This  flexibility  is  necessary  to  permit  the  axis  of  the  light 
beam  to  be  pointed  at  the  crack  on  the  test  specimen.  Image  distortion 
takes  place  if  the  crack  tip  is  located  near  the  edges  of  the  illumina¬ 
ted  region  of  the  specimen.  The  optical  platform  contains  a  15  mW 
continuous  wave  helium/neon  laser,  beam  steering  mirrors,  a  spatial 
filter  to  expand  the  beam,  a  lens  to  collimate  the  expanding  beam,  and 
finally,  a  second  lens  to  focus  the  beam  on  the  specimen  surface. 

The  loading  frame  is  powered  by  an  electric  motor-driven  screw 
jack.  It  features  a  rail-guided  loading  head,  and  a  series  of  springs 
to  permit  increased  accuracy  in  setting  the  desired  load.  Load  magni¬ 
tude  is  measured  by  a  series  of  interchangeable  calibrated  load  cells. 

The  magnitude  of  the  stress  intensity  factor  as  determined  by 
caustics  is  dependent  on  the  diameter  of  the  caustic  on  the  image  screen, 
and  on  the  optical  amplification  factor  chosen.  To  minimize  possible 
errors  due  to  these  two  factors,  the  following  procedure  was  adopted. 
First,  a  black  circle  of  known  diameter  was  applied  to  the  surface  of 
each  test  specimen.  When  illuminated  by  the  light  beam,  an  amplified 
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image  of  the  circle  is  seen  on  the  viewing  screen.  Since  for  these 
tests  all  measurements  were  t.o  be  made  from  photographs  of  the  caustic 
images,  the  recording  camera  scale  factor  was  obtained  by  photographing 
a  quality  steel  scale.  The  combination  of  recording  camera  scale 
factor,  and  the  ratio  of  the  circle  image  size  from  its  photograph  to 
the  actual  circle  size  on  the  specimen,  provides  the  optical  amplifica¬ 
tion  factor.  Representative  photographs  of  this  procedure  are  shown 
in  Fig.  3-3. 


0.160  INCH  DIAMETER  TEST  CIRCLE 


PHOTO  SCALE  FOR  SPECIMEN  7 


Fig.  3-3  Calibration  and  Photo  Scale  Shown  Above  Determine  the 
Amplification  Factor  >  . 


4.  DESCRIPTION  OF  EXPERIMENTS 


This  report  contains  the  results  of  testing  performed  with  seven 
1-in.  wide  test  specimens.  The  first  four  were  0.125-in.  thick 
2024-T851  aluminum  alloy.  The  fifth  specimen  was  0.250-in.  thick 
2024-T851  aluminum  alloy.  The  sixth  was  a  0.120-in.  thick  4340V  heat- 
treatable  steel  alloy  normalized  to  90  ksi  tensile  ultimate  stress. 
Finally,  number  seven  was  also  a  0.120-in.  thick  4340V  steel  heat 
treated  to  200  ksi  tensile  ultimate  stress. 

Specimen  geometry  is  defined  by  Grumman  Engineering  Drawing  number 
TGP-1075.  A  reduced  size  copy  of  this  drawing  is  shown  in  Fig.  4-1. 
Hole  faces  were  designated  A,  B  and  C  on  one  side  of  each  specimen,  and 
D,  E  and  F  on  the  other,  such  that  A/D,  B/E  and  C/F  were  faces  of  the 
corresponding  three  holes. 


Specimens  2-7  were  analyzed  for  both  the  cracked  open-hole  case, 
and  for  the  cracked  hole  with  an  interference  fit  fastener  installed. 
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DISCUSSION  OF  INDIVIDUAL  EXPERIMENTS 

Specimen  No.  1 

Material  -  2024  T851 

Width  -  1.00  in. 

Thickness  -  0.125  in. 

Hole  Dia  -  0. 188  in. 

Two  test  sequences  were  run  with  this  test  specimen.  The  first 
series  of  tests  was  conducted  on  the  specimen  in  which  the  central  hole 
had  been  cracked  to  0.090  in.  and  then  polished.  The  Kj  values  obtained 
by  caustics  were  50%  lower  than  that  predicted  by  the  theoretical  line. 
It  was  suspected  that  polishing  after  cracking  was  the  cause  of  this 
error. 

The  cracks  were  extended  for  all  three  holes.  The  crack  at  the 
central  hole  grew  from  0.090  in.  to  0.200  in.  into  previously  polished 
material.  Correlation  improved  somewhat.  Average  value  of  Ki  obtained 
by  caustics  was  77%  of  that  predicted.  Closer  examination  of  these 
first  caustics  photographs  indicated  that  the  apparent  caustic  diameter 
is  influenced  by  the  position  of  the  light  beam  relative  to  the  crack 
tip.  It  was  concluded  that  the  light  beam  must  be  centered  on  the 
crack  tip. 

A  re-run  of  the  center  hole  was  conducted  with  the  light  beam  cen¬ 
tered.  The  caustic  diameters  were  read  from  the  photographs  three  ways 
to  obtain  an  estimate  of  possible  scatter.  The  first  reading  employed 
the  Laboratory's  Fringe  Reader,  and  electro-optical  measuring  system. 

The  second  reading  was  obtained  by  smoothing  the  edges  of  the  caustic 
with  a  french  curve.  The  final  readings  were  obtained  by  smoothing  the 
edges  with  an  ellipse  template.  The  results  of  these  readings  are  pre¬ 
sented  in  Fig.  4-2. 

Finally,  the  upper  and  lower  holes  were  loaded  and  photographed. 
These  results  are  presented  in  Fig.  4-3. 
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STRESS  INTENSITY  FACTOR  -  KSI  x/lN. 


TEST  DATA  BAND 


(K|  -  ko  \/Li») 


K,  -  CAUSTICS  -  OPEN  HOLE 


SPECIMEN  DATA 


MATERIAL: 
WIDTH: 
THICKNESS: 
HOLE  DIA: 
CRACK  LENGTH: 
GEOM  FACTOR: 


2024-T851 
1.00  IN. 
0.125  IN. 
0.188  IN. 
0.201  IN. 
1.05 


22  24  26  28 
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GROSS  SECTION  STRESS  -  KSI 


R81-0244-002W 


REV  12-4-80 


Fig.  4-2  Stress  Intensity  Factors  (ksi  Caustics  vs  Theory,  Specimen  No.  1,  Center  Hole 


STRESS  INTENSITY  FACTOR  -  KSK  IN. 


GROSS  SECTION  STRESS  -  KSI 

R  8 1  -0244- 003 W _ 

Fig.  4-3  Stress  Intensity  Factors  (ksi  yin.), Caustics  vs  Theory, 
Specimen  No.  1,  Upper  and  Lower  Holes 
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Specimen  No.  2 

Material  -  2024-T851 

Width  -  1.000  in. 

Thickness  -  0.125  in. 

Hole  Dia  -  0.191  in. 

Open  hole  data  were  obtained  for  hole  faces,  A,  B,  E  and  F.  These 
four  crack  faces  had  crack  lengths  of  0.013  to  0.045  in.,  which  have 
geometry  factors  ranging  from  2.58  to  1.78.  The  Kj  values  obtained  by 
caustics  were  all  significantly  lower  than  predicted.  These  open  hole 
results  are  presented  in  Figs.  4-4  through  4-7. 

Steel  rivets  were  driven  into  the  three  holes  using  0.063-in.  thick 
aluminum  cover  plates.  Examination  of  the  crack  zones  indicated  that 
despite  the  use  of  cover  plates,  the  riveting  process  had  crushed  the 
surface  enough  to  mask  the  crack  and  the  crack  tip  caustics. 

The  cracks  were  extended  approximately  0.025  -  0.050  in.  Once 
again  the  specimen  was  examined  in  the  caustics  set-up,  but  the  caustics 
were  still  not  visible. 

The  cracks  were  extended  one  more  time  by  approximately  0.025  - 
0.050  in.  The  final  crack  lengths  ranged  from  0.008  -  0.162  in.  The 
caustics  produced,  where  visible,  were  of  a  quality  which  precluded 
their  use. 

Specimen  No.  3 

Material  -  2024-T851 

Width  -  1.00  in. 

Thickness  -  0.125  in. 

Hole  Dia  -  0.125  in. 

Open  hole  data  were  obtained  for  hole  faces  B,  C,  E  and  F,  which  had 
crack  lengths  of  0.134  -  0.168  in.  The  geometry  factors  ranged  from 
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STRESS  INTENSITY  FACTOR  -  KSI  s.m 


GROSS  SECTION  STRESS  -  KSI 


R81-0244-004W  REV  12-4-80 

Fig.  4-4  Strew  intensity  Factors  (ksi  Caustics  »s  Theory,  Specimen  2 A 
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STRESS  INTENSITY  FACTOR  -  KSI 


Fig.  4-5  Stress  Intensity  Factors  (ksi  VTrT.),  Caustics  ws  Theory,  Specimen  2B 
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STRESS  INTENSITY  FACTOR  -  KSI  V/IN- 


R81-0244-006W 


REV  12-5-80 


Fig.  4-6  Stress  Intensity  Factors  (ksi  >/■"•>.  Caustics  vs  Theory,  Specimen  2E 
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1.00  to  0.98.  As  shown  in  Figs.  4-8,  4-9,  and  4-10,  reasonable  correla¬ 
tion  between  the  predicted  values  of  Kj  and  the  values  obtained  by  the 
caustics  method  were  obtained  for  hole  faces  B,  C  and  E.  Hole  face  F 
was  very  different,  and  produced  a  Kj  curve  much  like  those  produced 
by  Specimen  2,  as  shown  in  Fig.  4-11. 
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Fig.  4-8  Stress  Intensity  Factors  (ksi  %/fn.).  Caustics  vs  Theory,  Specimen  3B 


STRESS  INTENSITY  FACTOR  -  KSI  v/lN. 
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SPECIMEN  DATA 


MATERIAL: 

2024-T8 

WIDTH: 

1.00  IN. 

THICKNESS: 

0.125  IN. 

HOLE  DIA: 

0.125  IN. 

CRACK  LENGTH: 

0.166  IN. 

GEOM  FACTOR: 

0.99 

RANGE  OF  POSSIBLE  VALUES 


_L 


-L 


10  12  14  16  18 

GROSS  SECTION  STRESS  -  KSI 


20 


22 


24 


R81-0244-009W 


REV  10-14-80 


Fig.  4-9  Strew  Intensity  Factors  (ksi  \ZlrTt,  Caustics  vs  Theory,  Specimen  3C 
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STRESS  INTENSITY  FACTOR  -  KSI  s'TN. 


GROSS  SECTION  STRESS  -  KSI 

R81  0244-010W _ REV  10-14-80 

Fig.  4-10  Stress  Intensity  Factors  (ksi  vln.),  Caustics  vs  Theory,  Specimen  3E 
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STRESS  INTENSITY  FACTOR  -  KSI 


CAUSTICS  -  OPEN  HOLE 


SPECIMEN  DATA 


MATERIAL: 

2024-T8 

WIDTH: 

1.00  IN. 

THICKNESS: 

0.125  IN. 

HOLE  DIA: 

0.125  IN. 

CRACK  LENGTH: 

0.168  IN. 

GEOM  FACTOR: 

0.98 

RANGE  OF  POSSIBLE  VALUES 


2  4  6 


8  10  12  14  16  18 

GROSS  SECTION  STRESS  -  KSI 


20 


22  24 


R81-0244-01 1W 


REV  10*14-80 


Fig.  4*11  Stress  Intensity  Factors  (ksi  >/Fn.),  Caustics  vs  Theory,  Specimen  3F' 
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Hole  faces  B/E  and  C/F  were  enlarged  to  0.191  in.  diameter,  and 
steel  rivets  with  an  initial  0.001  in.  diametral  interference  installed. 
The  rivet  in  hole  B/E  was  upset  using  the  6-in.  Stress  Wave  Riveter  at 
5.5  kilovolts  (kv).  The  rivet  in  hole  C/F  was  upset  at  7.0  kv.  The 
faces  of  the  specimen  were  protected  locally  with  1/8-in.  thick  rubber 
pads  during  the  riveting  operation.  With  the  rivets  in  place  large 
static  residual  caustics  were  observed  at  faces  B  and  F.  The  equivalent 
Kj  were: 

For  hole  face  B,  Kj  =  16,133  psi  /in. 

For  hole  face  F,  Kj  =  85,054  psi  /Tn. 

There  was  evidence  of  residual  caustics  at  hole  faces  C  and  E,  but 
the  patterns  were  not  sufficiently  clear  to  allow  the  Kj  values  to  be 
determined. 

The  rivets  were  removed,  and  it  was  noted  that  the  static  residual 
caustics  reduced  in  size  to: 

For  hole  face  B,  Kj  =  14,399  psi  /in. 

For  hole  face  F,  Kj  =  28,497  psi  /Tn. 

Finally,  hole  face  B  was  loaded  incrementally  to  a  gross  section 
stress  of  18  ksi,  and  the  Kj  values  plotted.  Figure  4-12  shows  the 
Kj  value  remained  constant  at  14,399  psi  /in’,  between  0  and  10  ksi,  then 
increased  linearly  to  23,157  psi  /Tn.  at  18  ksi. 

Specimen  No.  4 

Material  -  2024-T851 

Width  -  1.000  in. 

Thickness  -  0. 125  in. 

Hole  Dia  -  0.190  in. 

Open-hole  data  were  obtained  for  all  six  faces  on  this  specimen. 

The  crack  lengths  varied  from  0.039  -  0.060  in.,  with  geometry  factors 
of  1.85  -  1.58.  Above  4  ksi,  the  Kj  obtained  by  caustics  was  less 


r 


than,  or  in  one  case  (4D)  approximately  equal  to,  the  values  predicted 
theoretically.  These  results  are  plotted  in  Figs.  4-13  through  4-18. 

Steel  rivets  were  installed  in  holes  B/E  and  C/F.  The  rivet  in 
hole  B/E  was  driven  with  the  6-in.  Stress  Wave  Riveter  (SWR)  hit  on 
both  sides  at  5.5  kv.  The  residual  stress  field  surrounding  the  rivet 
was  high  enough  to  completely  mask  the  caustic  patterns.  The  rivet  in 
hole  C/F  was  driven  with  the  6-in.  SWR  hit  on  one  side  only  at  5.5  kv. 
It  was  observed  that  the  F  face  of  the  specimen  became  concave  (i.e., 
the  overall  field  reflection  became  smaller),  and  neither  the  crack  nor 


Fig.  4-12  Stress  Intensity  F actors  (ksi  \AirT.)  Caustics  vs  Theory, 
Specimen  3B  After  Removal  of  Rivet 
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STRESS  INTENSITY  FACTOR  -  KSIv'N. 


SPECIMEN  DATA 
MATERIAL:  2024-T851 

WIDTH:  1.00  IN. 


0  2  4  6  8  10  12  14  16  18  20  22 


GROSS  SECTION  STRESS  -  KSI 

R  8  1-0244-01  5W _ REV  10-4-80| 

Fig.  4-15  Stress  Intensity  Factors  (ksi  vTii-),  Caustics  vs  Theory,  Specimen  4C 
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STRESS  INTENSITY  FACTOR  -  KSI  VW- 


GROSS  SECTION  STRESS  -  KSI 

R81-0244-016W  REV  10-14-80 


Fig.  4-16  Stress  Intensity  Factors  Iksi  >/•”■).  Caustics  vs  Theory,  Specimen  4D 
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Fig.  4-17  Stress  Intensity  Factors  (ksi  \/TfO,  Caustics  ys  Theory,  Specimen  4E 
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STRESS  INTENSITY  FACTOR  -  KSI  v  IN 


Fig.  4-18  Stress  Intensity  Factors  (ksi  v/irT),  Caustics  vs  Theory,  Specimen  4F 


Specimen  No.  5 
Material  -  2024-T851 


Width  -  1.000  in. 

Thickness  -  0.250  in. 

Hole  Dia  -  0.191  in. 

Cracks  lengths  of  0.045  -  0.048  in.  were  obtained  on  one  side  of  the 
specimen  only  (i.e.,  these  were  corner  cracks).  Open  hole  caustics 
were  recorded  for  hole  faces  A  and  B.  Face  C  did  not  produce  a  caustic 
due  to  a  scratch  mark  on  the  specimen  made  during  handling.  For  a 
corner  crack,  the  theoretical  line  is  defined  by: 

Kj  =  ko  y Ltt 

where:  _ 

k  =  0.6557  /secant  jy(S  :x)] 

For  hole  face  A:  L  =  0.048  in.  D  =  0.190  in.  W  =  1.00  in.  k  =  0.6822 

For  hole  face  B:  L  =  0.045  in.  D  =  0.190  in.  W  =  1.00  in.  k  =  0.6813 

Figures  4-19  and  4-20  show  the  results  of  loading  this  specimen.  It 

is  shown  that  in  both  cases  the  slopes  of  the  Kj  curve  obtained  by 

caustics  were  less  than  the  slope  of  the  theoretical  line. 

An  aluminum  rivet  was  installed  into  the  A/D  hole.  The  riveting 
process  produced  an  epicycloid-shaped  zone  as  shown  in  Fig.  4-21.  The 
caustic  at  the  crack  tip  is  visible  along  the  cusp  axis.  The  caustic 
diameter  (0.45  in.  -  photo  size)  remained  constant  in  the  range  0-18  ksi. 

Specimen  No.  6 

Material  -  4340V  -  Normalized  to  90  KSI. 

Width  -  1.00  in. 

Thickness  -  0. 120  in. 

Hole  Dia  -  0.250  in. 
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STRESS  INTENSITY  FACTOR  -  KSI  s  IN 


A.  CAUSTIC  OF  RIVETED  AREA 


B.  EXPLANATORY  DIAGRAM  OF  CAUSTIC  SHOWN  ABOVE 
K81  0244-048W _ 

Fig.  4-21  Caustic  of  Riveted  Area  for  Specimen  No.  5,  Hole  A/D 
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This  specimen  was  the  first  of  two  steel  specimens  treated.  Crack 
lengths  of  0.028  -  0.047  in.  were  obtained  at  the  six  hole  faces. 

Loading  was  started  at  face  A,  and  it  was  noted  that  crack  tip  yielding 
occurred  between  28.  -  37.5  ksi.  The  Kj  by  caustics  was  higher  than 
the  theoretical  line  at  all  points  recorded.  The  data  available  is 
plotted  in  Fig.  4-22.  Large  residual  caustics  were  observed  at  the 
other  five  hole  faces.  The  Kj  range  was:  54.9  -  134.1  ksi^/in\ 

Steel  rivets  were  driven  into  all  three  holes.  The  heads  of  the 
rivets  on  faces  A  through  C  were  ground  away  to  improve  visibility. 
However,  surface  distortion  masked  all  crack  and  crack  tip  caustics. 

The  cracks  on  the  specimen  were  extended  by  additional  fatigue 
cycles  to  determine  if  the  caustics  would  reappear  as  follows: 

Hole  Letter:  A  B  C 

Initial  Crack  Length:  0.042  0.028  0.024 

Final  Crack  Length:  0.042  0.049  0.040 

Photographs  were  taken  of  the  entire  zone  between  the  edges  of  the 
specimen  and  the  rivet.  There  were  no  visible  crack  or  crack  tip 
caustics. 

Specimen  No.  7 

Material  -  4340V  -  Heat  treated  to  200  ksi 

Width  -  1.00  in. 

Thickness  -  0. 120  in. 

Hole  Dia  -  0.250  in. 

Crack  lengths  of  0.031  -  0.049  in.  were  generated  in  this  heat 
treated  steel  specimen.  The  resulting  geometry  factors  (k)  were  2.166 
to  1.878.  Loading  was  accomplished  over  a  range  of  7.5  to  37.5  ksi 
gross  section  stress.  In  all  five  readable  cases,  the  Kj  by  caustics 
was  substantially  lower  than  predicted  by  the  theoretical  line.  Figures 
4-23  through  4-27  show  these  results. 


GROSS  SECTION  STRESS  -  KSl 

R<Jl-0244-02lW 


Fig.  4-22  Stress  Intensity  Factors  (ksi  Caustics  vs  Theory,  Specimen  6A 
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STRESS  INTENSITY  FACTOR  -  KSI  v  IN. 


SPECIMEN  DATA 


MATERIAL.  4340V  -  200  KSI 


R81-0244-022W 


Fig.  4-23  Stress  Intensity  Factors  (ksi  VhT.),  Caustics  vs  Theory,  Specimen  7B 
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STRESS  INTENSITY  FACTOR  -  KSI 


Fig.  4-24  Stress  Intensity  Factors  (ksi  vTrT),  Caustics  vs  Theory,  Specimen  7C 
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STRESS  INTENSITY  FACTOR  -  KSI  VTO. 


Fig.  4-25  Stress  Intensity  Factors  (ksi  \/in.),  Caustics  vs  Theory,  Specimen  7D 
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SPECIMEN  DATA 


MATERIAL: 

4340V  -  200  KSI 

28  "  WIDTH: 

1.00  IN. 

THICKNESS: 

0.120  IN. 

HOLE  DIA: 

0.250  IN. 

26  -  CRACK  LENGTH: 

0.047  IN. 

GEOM  FACTOR: 

1.901 

24  - 

K.  -  ko  JGi. 

22  - 

1  N 

.  K,  -  CAUSTICS  -  RIVETED  HOLE 


- Q - B- 


0  4 

R8  1  -02  4  4-0  26 W 


r- 


.  k'.  —  CAUSTICS  -  OPEN  HOLE 


16  20  24  28  32 

GROSS  SECTION  STRESS  -  KSI 


40  44 


Fig.  4-27  Stress  Intensity  Factors  (ksi  Caustics  vs  Theory,  Specimen  7F 
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Steel  rivets  were  driven  into  all  three  holes.  The  cracks  and 
crack  tip  caustics  were  still  visible  in  all  five  holes  of  this  hard¬ 
ened  material.  For  gross  section  stresses  of  0-40  ksi,  the  caustic 
size  remained  constant  in  all  cases.  On  two  hole  faces,  B  and  C,  the 
loading  was  continued  to  70.8  ksi  (8500  lbs).  For  hole  face  B,  the 
caustic  diameter  increased  sharply  between  40-50  ksi,  then  remained 
constant  through  71  ksi.  Hole  face  C  showed  a  very  slow  irregular  in¬ 
crease  in  caustic  diameter.  With  the  exception  of  hole  face  D,  all  of 
the  holes  returned  to  their  initial  zero  load  diameters  after  loading 
to  71  ksi.  The  results  of  testing  on  Specimen  No.  7  are  plotted  in 
Figs.  4-23  through  4-27. 
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5.  DISCUSSION  OF  RESULTS 


The  purpose  of  this  report  was  to  carry  out  a  feasibility  demon¬ 
stration  of  setting  up  and  understanding  the  method  of  caustics  to 
determine  Kj  in  situations  where  the  results  of  the  conventional  ana¬ 
lytical  method  may  not  be  adequate.  We  believe  that  feasibility  is 
clearly  demonstrated  in  the  results.  Although  the  basic  relationship 
employed  in  relating  caustics  to  Kj  (Eq  (40)  of  Section  2)  requires 
the  value  of  D  to  be  measured  experimentally,  the  equation  itself  is 
based  on  the  theoretical  premises,  developed  by  Muskhel i shvi 1 i  and 
Westergaard  and  stated  by  Paris  and  Sih  (Ref  5)  and  therefore  the 
caustic  method  is  not  a  purely  experimental  one.  One  can  say  that  both 
equations  represent  a  model  of  stress  and  strain  at  the  tip  of  the 
crack,  with  the  caustics  equation  providing  a  more  detailed  representa¬ 
tion  of  actual  conditions  since  we  measure  the  local  conditions  in 
each  individual  case.  The  Paris  and  Sih  equation  (actually  first  de¬ 
rived  by  Bowie  (Ref  6):  Kj  =  ko/Gr,  implies,  for  example,  that  at  a  =  0, 
Kj  =  0.  Yet,  we  know  that  at  any  crack  tip  there  is  always  a  residual 
elastic  stress  field  that  assigns  Kj  a  value,  even  at  o  =  0.  The 
caustics  method  easily  infers  this  "zero  load"  condition. 

Other  difficult  conditions  exist  such  that  a  simple  theoretical 
equation  for  Kj  cannot  easily  be  derived.  One  such  case  is  the  exis¬ 
tence  of  residual  compressive  stress  around  the  hole.  Again,  in  this 
case  the  existence  or  lack  of  the  caustic  can  be  experimental  evidence 
of  the  effect  of  the  residual  stress  field. 

The  results  obtained  with  the  method  of  caustics  also  indicate  that 
the  opening-mode  stress  intensity  factor  Kj  is  not  a  linear  function  of 
load,  as  implied  by  the  equation  Kj  =  ka/Lrr.  We  do  not  have  an  ex¬ 
planation  for  this  as  yet,  but  we  do  believe,  as  stated  above,  that 


the  method  of  caustics  represents  the  behavior  of  the  crack  more 
accurately  than  the  classical  theory.  The  results  indicate  that  there 
is  greater  disagreement  between  the  method  of  caustics  and  the  equation 
Kj  =  ko/Lir,  for  1.3  <  k  <  2.00.  This  disagreement  occurs  with  small 
cracks  (i.e.,  0.005  in.  -  0.050  in.)  and  although  it  is  not  yet  proven, 
we  believe  that  it  is  due  to  the  biaxial  stress  field  near  the  hole, 
which  only  the  caustics  method  takes  into  account. 

The  tests  that  were  carried  out  with  the  steel  specimens  (6  and  7) 
showed  that  for  the  same  load,  the  caustics  of  the  normalized  (softer) 
specimen  were  considerably  larger  than  the  corresponding  caustics  of 
the  hardened  specimen.  This  behavior  was  expected  from  physical  rea¬ 
soning,  since  a  lower  yield  stress  will  allow  a  more  extensive  de¬ 
formation  by  plastic  deformation  prior  to  parting  at  the  crack  tip. 

Since  Eq  (40)  of  Section  2  was  derived  for  purely  elastic  cases,  it 
cannot  show  this  behavior,  nor  can  the  Paris  and  Si h  equation.  The 
sudden  increase  in  the  caustic  size,  however,  can  be  a  useful  indica¬ 
tion  of  incipient  plastic  deformation,  an  indication  not  easily  dis¬ 
played  by  analytical  models. 

With  the  above  principles  in  mind,  we  proceed  to  discuss  the  results 
presented  in  the  preceeding  section  in  more  detail.  In  several  of 
these  cases,  no  caustic  appeared.  Because  of  the  newness  of  the  tech¬ 
nique,  it  is  usually  not  possible  to  be  sure  what  the  reason  may  be 
for  lack  of  a  caustic  in  a  given  case.  It  may  indicate  that  Kj  has 
gone  to  zero  because  of  the  beneficial  effects  of  residual  stresses. 

It  may  be  that  the  caustic  has  been  destroyed  by  surface  distortions 
induced  by  the  riveting  process.  The  appearance  of  an  enhanced  caustic 
after  riveting  is  currently  interpreted  to  imply  that  the  crack  tip  is 
in  the  region  of  residual  tension  stress,  where  crack  propagation 
would  be  enhanced  by  the  residual  stress  field.  These  points  should  be 
kept  in  mind  in  the  following  discussion. 

In  total,  25  open  pre-cracked  hole  faces  were  analyzed  without 
rivets  installed  (open-hole  case).  A  review  of  the  plots  (Figs.  4-2 
through  4-27)  shows  that  in  only  one  of  these  cases  (6A)  was  the  Kj 


value  obtained  by  caustics  substantially  and  consistently  higher  than 
the  theoretical  line.  This  was  the  case  for  annealed  4340  steel,  and 
is  explained  by  excessive  plastic  deformation,  as  discussed  above.  In 
general  then,  the  test  data  tend  toward  the  conclusion  that  for  mean¬ 
ingful  fatigue  stress  levels  ( i . e . ,  for  gross  section  stresses  greater 
than  6-8  ksi)  the  K]  value  obtained  by  caustics  is  less  than  or  equal 
to  the  theoretical  line.  With  some  exceptions  (3F,  4B,  4C,  and  4D)  if 
geometry  factor  is  greater  than  1.3,  the  Kj  inferred  by  caustics  is 
substantially  less  than  the  theoretical  line.  For  geometry  factors  less 
than  1.3  the  Kj  obtained  by  caustics  is  approximately  equal  to  that 
obtained  from  the  theoretical  line. 

When  steel  interference-fit,  rivets  were  installed  in  the  aluminum 
or  soft  steel  specimens  significant  surface  distortion  occurred.  This 
distortion  is  a  natural  result  of  the  plastic  deformation  that  is  the 
goal  of  stress  wave  riveting.  In  only  three  cases  (3B,  3F,  and  4C) 
were  residual  caustics  visible.  These  cases  produced  large  static 
residual  Kj  values,  which  did  not  change  with  applied  load  through 
gross  section  stresses  equal  to  18  ksi.  The  crack  lengths  for  these 
cases  (0.134,  0.168,  and  0.060  in.,  respectively)  place  the  crack  tip 
outside  the  region  of  compressive  residual  stress  into  the  zone  of 
hoop  tension.  A  large  static  residual  Kj  would  be  the  expected  result, 
and  cracks  of  these  lengths  would  be  accelerated  by  riveting  if  our 
interpretation  is  correct.  The  absence  of  a  caustic  in  the  balance  of 
the  cases  implies  that  compressive  residual  stresses  work  to  reduce 
or  eliminate  stress  intensity  factors,  thus  working  to  arrest  crack  growth. 

When  steel  interference-fit  rivets  were  installed  into  the  hardened 
steel  specimen  (7B-7F)  static  residual  caustics  were  observed  to  exist 
at  all  five  readable  hole  faces.  These  caustics  remained  constant  in 
size  through  applied  loads  producing  gross  section  stresses  of  40  ksi. 

When  the  loading  was  continued  beyond  40  ksi  to  71  ksi,  the  caustics 
were  observed  to  grow  slowly  and  irregularly.  It  was  also  noted  that 
with  one  exception  (7D)  all  of  the  caustics  returned  to  their  initial 
zero  load  diameters  after  being  loaded  to  71  ksi  gross  section  stress. 

The  caustics  at  hole  "0"  enlarged  from  a  Kj  value  of  8.74  ksi  yiir.  to 
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16.18  ksi  in.  at  zero  load.  As  before,  it  remained  constant  through 
40  ksi,  then  increased  sharply  to  a  Kj  =  35.54  ksi  /Tr T.  at  70  ksi  gross 
section  stress.  This  interesting  behavior  is  deserving  of  additional 
study. 

It  should  be  pointed  out  that  clarity  of  the  caustic  epicycloid, 
and  the  ability  to  accurately  measure  it,  is  directly  related  to  the 
flatness  of  the  surface.  Due  to  the  exploratory  nature  of  this  first 
test  series,  all  of  the  specimens  were  polished  by  hand.  Better  polish¬ 
ing  methods  (i.e.,  machine  lapping)  should  be  used  for  any  subsequent 
work  in  this  field.  It  is  possible  that  some  or  all  of  the  anomalies 
seen  in  these  tests  are  the  result  of  distortion  of  the  caustic  epicy¬ 
cloid  by  localized  surface  irregularities. 

Specimen  7D  has  been  chosen  to  illustrate  the  caustics  method, 
since  it  represents  the  clearest  of  the  caustic  photographs  obtained, 
and  is  representative  of  the  test  procedure  which  was  followed.  Fig¬ 
ure  5-1  shows  the  caustics  for  the  cracked  open  hole  case.  Figure 
5-2  shows  the  same  hole  with  an  interference  fit  fastener  installed 
during  the  loading  series  in  which  all  the  holes  were  loaded  to  a 
maximum  of  40  ksi.  Finally,  Fig.  5-3  shows  the  same  hole  during  a 
later  series  of  loads  in  which  the  hole  had  already  been  loaded  to  a 
maximum  of  71  ksi.  These  three  pictures  are  of  caustics  for  which  Kj 
values  are  plotted  in  Fig.  4-25.  This  is  a  case  for  which  riveting 
had  little  beneficial  effect  in  reducing  Kj,  probably  because  the  high 
yield  stress  prevented  significant  plastic  deformation  in  the  surround¬ 
ing  steel,  so  the  beneficial  residual  compressive  stress  did  not  have 
much  effect.  In  those  cases  where  the  benefits  of  residual  stress  were 
great,  the  caustics  are  either  non-existent  or  much  less  detectable 
than  in  specimen  7D. 
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Fig.  5  3  Specimen  7D  -  with  Rivet  -  Later  Loading  to  70  ksi  (Sheet  1  of  2) 
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6.  CONCLUSIONS 


The  method  of  caustics  promises  to  be  a  very  useful  tool  for  the 
investigation  of  stress  intensities  in  local ized-stress  regions  such  as 
those  near  crack  tips.  As  a  complement  to  purely  analytical  treatment 
of  such  regions,  it  provides  direct  measurements  of  surface  distortions 
in  such  regions,  giving  insight  into  such  factors  as  the  influence  of 
imposed  residual  stress,  the  natural  residual  stress  field  induced  as 
the  crack  forms,  the  effects  of  hardness,  and  perhaps  other  factors  not 
yet  discovered. 

This  preliminary  investigation  of  the  use  of  caustics  indicates 
that  stress  intensity  factors  follow  the  general  characteristics  of 
theoretical  predictions,  except  in  certain  special  situations.  At  low 
stress  levels,  the  caustics  results  indicate  a  zero-load  value  of  Kj 
that  is  substantial,  and  not  predicted  by  the  classical  (elastic)  theory. 
At  values  of  the  geometric  factor  greater  than  1.3,  the  caustics  values 
are  substantially  below  the  theoretical  predictions.  As  loading  begins 
to  approach  the  yield  stress  of  the  material,  there  is  a  tendency  for 
the  caustic  values  to  increase  relative  to  the  theoretical  predictions. 

When  plastic  deformation  of  the  region  surrounding  the  precracked 
hole  is  introduced  by  installation  of  a  stress  wave  driven  rivet,  the 
caustics  disappear  if  the  tip  of  the  crack  lies  within  the  region  of 
residual  compressive  stress;  at  least  until  the  loading  is  increased  to 
high  values.  In  a  few  cases  when  the  crack  tip  lies  outside  the  com¬ 
pressive  residual  region,  the  no-load  Kj  indicated  by  caustics  is  sub¬ 
stantially  increased.  Preliminary  indications  from  these  results  are 
that  stress  wave  installed  rivets  can  arrest  cracks  up  to  a  critical 
length  that  depends  on  the  rivet  installation  parameters,  at  least  for 
loading  less  than  a  particular  level. 
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FRBC  SOI  N3  F AOS  BUNK-MOT  II 


APPENDIX 


FRACTURE  MECHANICS  BACKGROUND 

Structural  materials  of  crystalline  structure  usually  fracture  at 
strengths  well  below  their  theoretical  strength  levels.  Such  premature  fail¬ 
ures  are  caused  by  the  many  structural  imperfections  which,  when  combined 
with  oscillating  tension  loads,  initiate  cracks  that  result  in  fractures. 

The  contemporary  approach  to  structural  design  accepts  material  imper¬ 
fections  as  facts  of  life.  To  accommodate  their  presence  and  avoid  their 
catastrophic  effects,  designers  and  stress  analysts  use  various  theoretical 
and  empirical  techniques. 

Fracture  analysis  identifies  three  types  of  cracks  as  shown  below: 


Our  applications  relate  to  the  cracks  emanating  from  a  hole,  or  Mode  I  crack 
surface  displacement,  as  shown  in  Fig.  A-l. 
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WHERE  S  =  REMOTELY  APPLIED  TENSION. 

Fig.  A-1  Specimen  Containing  Mode  I  Crack 


where  S  =  remotely  applied  tension. 

Dimensional  analysis  and  geometrical  considerations  show  that  the  stress 
intensity  factor  for  Mode  I  is  given  by 

KI  = 

Following  Ref.  5,  the  elastic  stress  distribution  at  the  crack  tip  for 
Mode  I  is  given  by  the  following  equations  (see  Fig.  A-1.) 


0  r  .  0  .  30 1 

—  1  -  sin-  sin-r- 


Txz  Tyz 


and  ~ 


!=  0  pla 

V  (ox  + 


ane  stress 


o  )  plane  strain 

y 


where  o  ,  o  ,  a  are  stresses  in  the  x,  y,  Z  directions  and  t  ,  x  , 
x  y  z  xy  x  z 

x  are  shear  stresses. 

and  Kj  =  the  stress  intensity  factor. 

The  equations  above  clearly  show  that  for  any  radius  r,  the  intensity 
of  the  stresses  depends  on  the  value  of  Kj.  For  example,  taking  an  element 
at  Q  =  0  will  result  in  o  =  a  =  KT/  J  2irr. 

The  crack  distribution  is  as  shown  below. 


The  meaning  of  is  closely  related  to  that  of  K^,  the  stress  concentra¬ 
tion  factor.  Kj  and  Ky,  however,  are  not  interchangeable.  For  example,  the 
hoop  stress  distribution  of  a  uniaxially  loaded  specimen  with  a  hole  in  its 
center  is  shown  below: 
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In  crack  problems  with  the  crack  embedded  in  an  infinite  plate  at  an 
angle  3  with  the  axis  of  loading,  we  can  use  the  complex  stress  function 
v  (c)  for  the  region  close  to  the  crack  tip.  From  Ref  5,  the  stress 
function  is  given  by 


(c) 


K* 

2(2ttc)1/2 


(A- 


where  K*  =  Kj  -  iKjj 

Kj  =  stress  intensity  factor  for  the  opening  mode 


Ktt  =  Stress  intensity  factor  for  the  sliding  mode 


